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The method of continuous canonical transformation is apphed to the double 
exchange model with a purpose to eliminate the interaction term responsible 
for non conservation of magnon number. Set of differential equations for the 



o 

Q ' effective Hamiltonian parameters is derived. Within the lowest order (approx- 



imate) solution we reproduce results of the standard (single step) canonical 
transformation. Results of the self consistent numerical treatment are com- 
pared with the other known studies for this model. 
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I. INTRODUCTION 

The discovery of colossal magnetoresistance in the doped inanganites attracted new inter- 
ests in studying itinerant ferromagnetism phenomena. The double exchange model (DEX), 
introduced a long time ago [0,0, seems to be a good starting point to explain the param- 
agnetic - ferromagnetic transition. Conduction Cg electrons interact via the strong Hund's 
coupling with the localized Mn ions (their spin being S = 3/2). This interaction drives the 
core Mn spins to ferromagnetic alignment, owing to the kinetic processes of itinerant elec- 
trons. Whether the DEX scenario itself is enough for an explanation of magnetoresistance in 
manganites, or should be supplemented by other realistic effects like e.g. lattice Jahn- Teller 
distortions P], competition with the superexchange 0,^ and strong Coulomb interactions 
between electrons 0, is an open question under discussion. 

In this paper we want to reexamine physics of the DEX model using the renormalization 
group approach in a version proposed recently by Wegner and independently by Glazek 
and Wilson j^. This procedure known as the flow equation method has proved to be a 
powerful tool, especially for analysis of models composed of several coupled subsystems. So 
far, it has been successfully applied to the problem of electron - phonon interaction (where 
the Frohlich transformation has been revised), the single impurity Anderson model |T0[, the 
spin - boson model for dissipative systems |]TT|, the spin - polaron coupling for t—J model []T^|, 



DEX model in the RKKY (small coupling) limit JTSl, and the charge exchange interaction 



for the boson - fermion model |T^. The same method has also been used for studying 



effects of correlations, e.g. in the Hubbard model fTSl {t/U expansion), the large spin 
Heisenberg Hamiltonian [|l^ {^/S expansion), etc. Recently there has been proposed a highly 
sophisticated computer aided perturbation method based on the flow equation technique to 
study the dimerized spin models [|17|. It is our belief, that the flow equation technique can 



be a reliable source of information also for the DEX model. In particular we would like to 
consider the strong (ferromagnetic) coupling limit (relevant for manganites) and compare 
our results to other studies of this model. 



The DEX model for ferroinagnetism is described by the Kondo type Hamiltonian |T^ 

H = Y,{ek- f^) cl^Cka - Jh Y. (SiSi)'""' cLcia' (1) 

k,cr i,(T,cf' 

where c^l correspond to annihilation (creation) operators of the conduction Cg electrons, 
Sj denotes their spin operators and Sj stands for the spin operator of Mn ions. Local 
ferromagnetic interaction is characterized by the Hund's coupling Jh > 0, for manganites 
known to be very large. With a help of the Holstein-Primakoff transformation we can 
represent the spin operators Sj via the magnon operators a|, Oj such that S^^ = ai{2S — 
ajaj)-'^/^, S^ = {Slj^y, 5/ = S' — ajoj. Magnon operators a|, Oj obey the boson (commutation) 
relations. 

At sufficiently low temperatures the system is close to ferromagnetic (ground state) 
ordering, so we can simplify spin operators to S~ ~ al\^2S, S^ ^ aiV2S. Using the Pauli 
operators for Cg electron spins one gets the model Hamiltonian in a following form 

H = E^kcLcka + ^ E ai+q^P (4-qTCkT - 4-qiCki 
k,CT q,P,k 



-^^l/^EK^-qTCki+h-C.) , (2) 



S 

q,k 

with ^k = ^k — fJ'" and fi^ = ft + \SJh, fi^ = ft — \SJh- If temperature is comparable 
with Jn/kB then one should include the higher order terms in the square root expansion 



1 — {a\ai/2S) for magnon operators. 
In what follows, we design the unitary transformation to eliminate last part of the Hamil- 
tonian (0) linear in a'^P operators, which is responsible for a violation of the magnon number 
Y^i < ciidi >■ To the leading order of 1/5* one can eliminate this exchange interaction using 
a single step canonical transformation e^He~^ with the generating operator A given by |T^ 



However, this transformation generates a whole lot of higher order interactions, their ampli- 
tudes being eventually not negligible. Recently it has been shown by means of the standard 



perturbation treatment ||20| that these interactions give rise to quantum corrections for spin 



wave spectrum (both for dispersion and the hfe time effects). That such corrections are in 
fact important it is known independently from the direct numerical studies of finite chains 



21| and from analysis based on the variational wave functions p2 |. 

Instead of the single step transformation (^) we propose in this paper a different method 
using an infinite sequence of the infinitesimal transformations what gives us more control 
for a derivation of the required effective Hamiltonian. In particular, we want the higher 
order many-body interactions to be as small as possible, thus being more tractable via the 
standard perturbation study. 

In the following section we give a brief introduction to the method of continuous trans- 
formation and derive the corresponding flow equations for parameters of the DEX model. 
Next, we discuss an analytical approximate solution of these equations and compare it with 
results of the standard single step transformation. In the last part we present the selfconsis- 
tent numerical solution for the model parameters along with some rough estimation of the 
spin stiffness coefficient. 

II. FORMULATION OF FLOW EQUATIONS 

A main idea of the method is to transform the initial Hamiltonian H through the series of 
unitary transformations H{1) = U{l)HU'^{l), labeled with some continuous flow parameter 
I. In a course of transformation the Hamiltonian evolves according to the following flow 
equation 

^ = [V{1),H{1)], (4) 

where the generator ri{l) is related to U{1) via ri{l) = {dU{l)/dl)W{l). This operator has 
to be chosen depending on a purpose of the transformations. Wegner has shown ||^ that 
with ri{l) = [H{1) — Hint{l) , H {[)] one can eventually eliminate the (perturbation) part of 
the Hamiltonian Hintil -^ oo) = 0, provided that no degenerate states are encountered. 



Alternative choice for rj, efficient even in a presence of degeneracies, has been proposed 
recently by Mielke |2^. In this work we use the slightly modified Wegner's proposal for ri{l). 
To be specific, we define the interaction part as 

Hintil) = ^ E (^k,q(/)a;cLq,TCki + h.c.) . (5) 

The remaining part Hq{1) = H(l) — Hint{l) may contain not only the other two terms of (0) 
but additionally also contributions induced by transformation for / > 0. To take these into 
account we assume Hq{1) to have the following structure 

^0(0 = E^k(OcLcka + ^ E ^k+q,k'+q' [f/k,q,q',k' (OCkiCqtCq/^Ck'i 
k,cr k,k',q,q' 

+ (^k,k',q,q'(04TCk'T - Mik',q,q'(04iCk'i) a\a^] + dH^) , (6) 

where 6Hq{1) contains all types of interactions not shown explicitly in (P). The initial 
conditions for the model parameters read 

e^(0)=£k-/i% f/k,q,q',k'(0) = , (7) 

WO) = J^^|, MkVq,q.(0) = ^ . (8) 

We choose the generating operator ri{l) = [J2k,a ^klOcko-Cko-, HintiJ)] which explicitly is given 
by 

-1 

k,q 



^(0 = -J^ E "k,q(0 (4,q(0«i^k-q,TCki - h.C.) , (9) 



where Q;k,q(0 = 'Ck-q(0 ~ ^k(0- Notice, that (^ has similar structure to the generating 
operator A of the standard transformation (|^). 
Using the general flow equation (^ we obtain 

dH{l) _ 1 ^^^ /,^^2^r /n.fJ .. , , i. . ^ A ^^ _ (n\T (n\2j 



Y. ("k,q(/)) Uk,q(0ai4-q,TCki + h.C.) " ^ ^ "k,q(0 |4,q(0 P^jCki 



V ^ k,q k,q 



^ k,k',q,q' 



^ E ("k,q(0 + ak',q'(0) 4,q(0^k',q'(0 [4^^' (4-qTCk'-q'T'^k,k' 
k,k',q,q' 

~ %J.CkJ.5k_q,k'-q' j + Ck'iCi^-qt^k'-q'TCki^q.q' 

+ [77(0, 5H,{1)] + 0(/M) + 0(/f/) . (10) 



Terms of the order I{l)M{l), I{l)U{l) are symbolically denoted as 0{IM) and 0{IU). If 
they were included in the diagonal part (||) through 6Hq{1) they would induce some higher 
order interactions given by [ri{l),6Ho{l)]. 



Equation ([1^) is a differential flow equation for the Hamiltonian which has to be solved. 
In the next section we solve it approximately neglecting the last three terms on the right 
hand side. 

III. LOWEST ORDER SOLUTION 

It is instructive to study first the flow equation ([T0|) with the last three terms on the right 
hand side omitted. It means that we neglect the interactions expressed by more than four 
operators, i.e. scattering between more than two particles. On a level of this assumption 
one may expect that effective Hamiltonian should differ from its initial form (^ by some 
correction comparable with H2 of the Ref. [^. The set of the flow equations for parameters 
of the DEX model Hamiltonian is simply given by 

^^'"^^^ - -«L(0/.,.(0 , (11) 

-lE«k,.(/)i/M(or, (12) 

q 

(ak',k'-p(0 + "k,k-p'(0) 4',k'-p(0-^k,k-p'(0 5 (13) 

("k+q,q(0 + ak'+q',q'(0) 4+q,q(0^k'+q',q' (0 > (1^) 

^l («k',q(0+«k,q'(/))/k',q(/)/k,q'(/)- (15) 

There is no renormalization for spin a =1 electrons, so ^k(/) = ^k. 

Still, the set of flow equations (|ri| -|T5|) is impossible to solve by other than numerical way. 
To get some insight in a process of the continuous transformation for H{1) we assume further 
that the energy ^k(/) is only weakly affected by renormalization. We verified validity of such 
assumption by solving selfconsistently the flow equations numerically for one dimensional 
tight binding electron dispersion (see the next section). We thus can drop / dependence 
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dl 


^iii) 


dl 


C^f^k.p.p'.k' (0 


dl 


<k',q,q'(0 


dl 


<k',q,q'(0 



of the parameter ak,q(0 '^^ ^^e right hand side of flow equations. One easily obtains an 
analytical solution for a flow of the exchange coupling 

/k,q(0 = /k,q(0)e-<.' = J^^e-(«-.-«>^)'' . (16) 

In the limit / -^ oo the exchange coupling drops asymptotically to zero and so does the part 
(|) of the Hamiltonian, Hintil — > C)o) = 0. 

Determination of the other /-dependent parameters, such as ^^{1), U{1) and M°"(/) is 
straightforward. We summarize our results by showing values of these parameters in the 
limit / -^ oo 

e,i(oo) = ei-4^E7r^, (17) 

^k,p,p',k'(oo) = — — /p,k',p',k (18) 

^k,k',q,q'(00) = Y" + %^ /k,k+q,k',k'+q' , (19) 

Mik',q,q' (00) = ^ + ^ /k'-q,k',k-q',k , (20) 



where 



f — ^1 '^2 + ^3 Ja (r,-,\ 



Using the standard canonical transformation (|^) one obtains the same scaling for the pa- 
rameters (|17|-|^) but with a different factor / [|20| 



/12m - \ (^ + ^) . (22) 

A general feature of the continuous canonical transformation is that it derives the effective 
Hamiltonians avoiding any singularities for the renormalized energies and interactions [0,§. 
For instance, the effective retarded interaction between electrons in a coupled electron- 



phonon system has been shown to be —\M^ ^^P^jq (^k+q — ^v) + 



< 



-1 



instead of the 



-1 



divergent Frohhch result I Mq' ^^pcUq (^k+q — ^k) ~ ^q W^ i^f ^'^ denotes the electron 
phonon coupling and cUq, ^k refer to phonon and electron energies respectively). 



A similar situation takes place in the DEX model studied here. For small values of JhS 
(as compared to the electron bandwidth) we obtain less divergent factor ( pT| ) than it has 
been predicted from the standard canonical transformation (^) found in the Refs |T9| , pO 



For the hmit of large JhS both factors (|2T| , |22| ) asymptotically approach / ~ —{JhS) ^ and 

effectively M'^ ~ 0, f/ Jj//2. 

In the limit / ^ oo, the exchange interaction @ is absent and then one can roughly 
estimate the magnon dispersion as 

^q = ^ E Kk,q,q4 - Ml^c^^^ni] . (23) 

Here the expectation value n^ = (c]^^C]ia-) has a meaning of the Fermi-Dirac distribution 
function of the argument (,^{1 = oo). By inspecting the expressions (plip2l) one notices that 
the diagonal terms /i, 2,3=1,4=2 are in both cases identical, independently on magnitudes of 
Jh and S. Effectively, the magnon spectrum becomes 

^q = ^E K -^k) + ^^E ,T A ■ (24) 

^^^ k ^^^ k ^k~Sk+q 

Our lowest order estimation ( p4D is thus in agreement with predictions based on the standard 
canonical transformation ||19|,^ and other earlier studies of this model [p5|-|27[] as well. The 
gapless Goldstone mode c^q=o = (for arbitrary temperature and Jh) marks the spontaneous 
breaking of rotational symmetry. 

IV. SELFCONSISTENCY CORRECTIONS 

In this section we take into account the effect of the terms neglected by us so far in the 
equation (|T0|) . We postulate that the effective Hamiltonian should be given in the form Hq{1) 
(^) with some small correction expressed there by 6Ho{l). Following other studies based on 
the flow equation method we consider effect of the higher order interactions by reducing 
them to normally ordered form. As can be seen below, 6Hq{1) would then be expressed by 
the fluctuations around the mean field values, which should be small. 
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i) Correction of the order 0{IU) arises from the commutator between 77 (^ and the 
electron electron interaction. It is easy to verify that this term can be expressed as follows 
1 



0{IU) 



Ar /ArH"k,q(0-^k,q(0 X! ^k',p',p",k"(0 '^k'+p',k"+p" 

rsyrs i^q k',p',p",k" 

X a\ (cp/|4/|Ck"iCki(5k-q,p" - 4-qTCp'tCp"TCk"i^k,k') + h.c. 

'■ 0{IU) : + T^7T^H^cA-^^^^iYl ["k',q(04',q(/)f/k',k-q,k'-q,k(0 
k,q k' 



A^VA^ 



^k'-q 



-^y - ak,q(04,q(0 (f^k,k'-q,k'-q,k(0^k'-q " ^k',k-q,k-q,k' (O'^k' ) + ^-C (25) 

ii) The other term 0{IM) comes from the commutator between rj and the magnon 
electron interaction. Its contribution to the flow equation (^ is 

0{IM) = -i= 5: ak,q(0^k,q(/) E [K,n" (^k',k",q',q"(04'TCk"T 

I\y/I\ k,q k',k",q',q" 

- ^k',k",q',q"(OCk'iCk"i) a^'Cl^-qtCki + a^ai'V (^k',k",q',q"(0Ck'TCki4-q,k" 
+ ^k',k",q',q"(0Ck-qTCk"i'5k,k' 



<^k'+q',k"+q" 



h.c. 



E <,k',q,q(04' 



: OilM) : + — ^ E «W-cTCki «k,q(/)/k,q(0 

NVN k,q I 

-^k',k',q,q(04') + E ^q' (A'^k-q,k-q,q',q'(0 + ^^k.q'.q'i 
q' 

+ E [(l - 4-q' + ^q') ^k-q,k-q',q,q'(0«k,q'(04,q'(0 
q' 

+ (4-q' + ^q') ^k-q',k,q,q-q'(0"k-q',q-q'(04-q',q-q'(0] } + ^.C. 



(26) 



'^q'^q/' 



We introduced here the expectation value for the magnon's number operator n^ 

which can by approximately taken as the Bose-Einstein distribution function of the magnon 

energy cUq. 

From now on, we put the both contributions ( p5| , |26| ) into the 6Hq{1) part of the Hamilto- 
nian and neglect there the normally ordered parts : 0{IU) :, : 0{IM) : which are supposed 
to be small. We neglect also a contribution to the flow equation which might arise from the 
commutator [rj, SHq]. This is the only simplification we need in order to close the set of flow 
equations for the Hamiltonian Hq{1) + H-intil). By inspecting the structure of expressions 
(^,^) we conclude that the flow equations (|12|-^) remain unchanged. The only quantity 



affected directly is tlie excliange coupling constant /k,q(0- The revised flow equation ([TT| ) is 
now given by 



^k' 



(. k' 

- (t/k',k-q,k-q,k'(/) - Mi[,^i,,_q^q(/)) 4,] - ^E^q' (^k-q,k-q,q',q' (0 

q' 

+ Mi[^k,q',q'(0)} + T7E«k',q(/)/k',q(0f^k',k-q,k'-q,k(0 (4'-q " 4') 
^^ k' 

- 17 E [(1 - 4-q' + ^q') ^^k-q,k-q',q,q'(0«k,q'(04,q'(0 

q' 

+ W-q' + ^q') ^k-q',k,q,q-q'(0«k-q',q-q'(04-q',q-q'(/)l • (27) 



We studied numerically the system of coupled flow equations (|T2]-^,^) solving them self- 
consistently via the Runge Kutta algorithm. Since the model parameters such as f/k,q,p,k+q-p 
and ^kqpk-q+p depend on three momenta it is a rather cumbersome task to study the 3 
or even 2 dimensional systems. For a demonstration we therefore used the one dimensional 
tight binding lattice with its initial dispersion 6\f_{l = 0) = —2t cos ka. From now on, we set 
the lattice constant a = 1 and choose the initial bandwidth as a unit W = 4t = 1 (flow 
parameter / will be expressed in units of W'"^). 

We discretized the first Brillouin zone with a mesh of 200 equally distant points. Starting 
from the initial conditions (0,|]) we computed iteratively the renormalized model parameters 
using a following scheme y{l + 61) = y{l) + y'{l)Sl where y'{l) = dy{l)/dl and it is given by 
one of the flow equations (|12|-|1^, ^7]) for a corresponding parameter y. Increment of the flow 
parameter was taken 61 = 0.0001 for / < 0.1 and 61 = 0.001 for / G (0.1; 1.0). 

Figure 1 shows how the model parameters evolve with an increasing /. Practically, 
already from / = 0.3 these quantities start to saturate at their asymptotic values. The 
exchange interaction disappears very fast, whereas the magnon - electron interaction is 
reduced roughly 10 times. There occurs some renormalization of ej^ electron energies and a 
simultaneous induction of a relatively strong electron - electron attractive interactions. 
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FIGURES 




u 



0.4 



/ 



FIG. 1. Evolution of the DEX model parameters with respect to a varying flow param- 
eter / at T = and hole concentration 0.3. Individual lines correspond to the follow- 
ing quantities: / = j^t Ek,q l-^k,qP, e = ^Ek(ek^)^ U = j^b- Ek,p,q |f^k,p,q,k+p-qP and 



^ = W Sk,p,q l^k,p,q,k-p+ql ' 




FIG. 2. Plot of the normalized exchange coupling /kq(^ = l)/-^kq(^ = 0). This interaction is 
weak enough to be considered as negligible. 
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Figure 2 illustrates what is left of the exchange interaction at / = 1. The highest values 
of /k,q are 7 x 10^^ of the initial interaction J^J S/2. Such small magnitude is in our 



opinion negligible and thereof we represent below the effective model parameters (formally 
corresponding to / = oo) through their values obtained at / = 1. 

Let us next have a brief look on the electron part. In figure 3 we plot the dispersion 
for spin a =| electrons obtained from the selfconsistent solution of the flow equation (|12|). 
Effective bandwidth becomes of the order ~ 0.83 of the initial bare bandwidth W . This 



renormalization is weaker than a result of the standard canonical transformation |119| , |20 
which gave the scaling factor ~ 0.64. As concerns the effective mass we notice a similar 
tendency. The whole band of a =| electrons drifts away from the partly occupied band of 
a =1 electrons increasing the gap between them which initially was JhS. The flow equation 
method gives a somewhat smaller shift than the standard canonical transformation. 




0.0 - 



-0.5 



cont. transf. 



free 



-71 ° k '^ 

FIG. 3. The effective dispersion for spin a =[ electrons obtained from a selfconsistent solution 
of the flow equations (solid line) and via the standard single step transformation (dashed line). 
The dotted curve shows the initial bare dispersion ek = —2t cos ka. 

Figure 4 shows a strength of the induced electron - electron interaction for the BCS 
channel f/k,-k,-q,q (scattering between the electron pairs of the total zero momentum). This 
interaction is again weaker by almost 30 % from the corresponding magnitude determined by 
the standard transformation. It should be stressed here that electron - electron interactions 

12 



may eventually play important role in a low energy physics of the DEX model only when hole 
concentration 1 — ra is close to zero. Otherwise, the scattering processes between particles 
originating from so much separated (~ JhS) electron bands would not become efficient. 



'-^k.-K-q.q 




FIG. 4. The effective interaction between electrons in the zero momentum BCS channel, i.e. 
f^k,-k,-q,qcLc_j^|C_q|Cqj. Other elements of the potential C^k,p,q,k+p-q take on average the same 
values as shown in this picture (around —0.8). 



If the initial Hunds coupling Jh is large and hole doping is not close to zero then the 
low energy physics of the effective Hamiltonian H{1 = oo) is mainly determined from the 
electron - magnon interaction. This interaction is characterized by two coupling constants 
^kpqk-p+q which, accordiug to our numerical estimation, can vary between —0.5 and 0.5 
depending on the momenta q, k, p and only sligthly on temperature T and concentration 
n. In general, in this method we find that absolute values of both coupling constants M'^ 
become reduced by 10 to 20% as compared to the result of the standard transformation 
(given in equation (1T9|) with the Z*^*^^ factor (p2l)) . Figure |^ shows this effect on the example 
of the diagonal part M^^^ ^ (notice, that in the approximate solution discussed in section 
III there was no difference between both methods for such diagonal elements). 
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Mk,M,q 



FIG. 5. Potential of the electron - magnon interaction Mj[ i^ „ „ obtained from a selfconsistent 
solution of the flow equations (solid lines) and from the standard canonical transformation (dashed 
curves). We used three representative values for the momentum q as marked on a right corner of 
the figure. 

A natural expectation is that reduction of the electron - magnon interaction would in a 
consequence affect such quantities like magnon spectrum, the life-time of these quasiparticles 
and finally also the Curie critical temperature. It is not a scope of this paper to study 
all these effects (we mainly intend to present this new technique for a derivation of the 
effective Hamiltonian). To get some insight we present below a rough estimation of the 
magnon dispersion based on the lowest order formula (pSf ). Higher order corrections are 
straightforward to carry out but we believe they would not alter anything in this context. 

In the lowest order perturbative determination of the magnon dispersion cuq we simply 
need only the diagonal parts M^y_^ ^ and at low temperatures main contribution comes from 
cr =^. In the long wavelength limit q ^ for cjq one gets a parabolic momentum dependence 
with the spin stiffness coefficient Dg = lim.q^o uj^/ q"^ . Using the standard expression ( |2^ 
Furukawa [^ gave an explicit formula for Dg at low T and studied Dg with respect to Jh 



and the hole concentration p = 1 — n (where n = I]k,o- < cj^o-Cko- >)• For a comparison we 
show in the figures || and |^ our results. 



14 




t/J 



H 



FIG. 6. Spin stiffness Ds = 2SDs/t as a function of the inverse Hund's coupling for T = and 
hole concentration 0.3. Solid line shows a result obtained from the flow equation method and the 
dashed one refers a standard estimation based on equation (M). 



D. 



0.01 - 



0.005 



Standard 
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FIG. 7. Spin stiffness D^ as a function of hole concentration for J///VF = 2, T = obtained 
from the flow equation method (solid curve) and the single step transformation (dashed line). 



15 



Due to a discussed above reduction of the electron - magnon coupling (see figure ||) we 
notice a partial softening of the spin stiffness in a whole regime of carrier concentration. 
Such a softening takes place mostly for the strong Hunds coupling. Going towards the limit 
of small Jh this effect becomes less pronounced. In particular, for Jh = '2W and hole 
concentration p = 0.3 our result for Dg is almost 20% smaller as compared to the standard 



result reported by Furukawa |26 



This tendency for a softening of the spin stiffness (especially in a limit of large Jh) agrees 
qualitatively with the results reported recently by Shannon and Chubukov @] who introduced 
a novel large 5* expansion scheme for systems with strong Hund's coupling. Authors showed 
that quantum effects caused a relative softening of spin-wave modes at the zone center. 
There is a rich literature where authors report differences between the Heisenberg cosine 
dispersion (for Uq) and the actual dispersion found from the calculations for DEX model 
P^ , |26|p8| . In those papers a relative flattening of the dispersion cjq near the zone boundaries 
has been found. Our results follow closely the same behavior. 

At the end of this section we would like to make an effort to analyze in more detail the 
effect of some neglected terms : 0{IU) : and : 0{I M) :. As an illustration, let us consider 
one of possible contributions to SHq{1) 

SH^^\l) = ^-jf= Y. ^q,k,p',p"(0 : aqCk-qTcl,'TV'TCk+p'-p"i : (28) 

NVN q,k,p',p" 

taken from : 0{IU) : (|25|) . Its initial (/ = 0) amplitude is of course Vq,k,p',p"(^ = 0) = 0. 
The flow equation corresponding to this potential is given by 

~Tf = — «k,q(0-^k,q(Of^k,p',p",k+p'-p"(0 • (29) 

We checked numerically that the potential V of the interaction (^Sf ) is negative. On aver- 
age, its value is 1/A^^ I]k,q,p',p" Ki,k,p',p" — —0.4 which is close to strength of the electron - 
magnon interactions < M"" >. To be speciflc we show in flgure 8 this potential in two chan- 
nels: Vo,k,-k,q(04ckTcikTC-qTCqi (top plcture) and V;j,k -k -k(0'^qCk-qTcikTC-kTCki (bottom 
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picture). 



V 



0,k,-k,-q 




'k,q,-q,-q 




FIG. 8. Potential of the many body interaction given in the equation (|2^) at T = and hole 
concentration 0.3 

Such many body interactions would however be not efficient at low temperatures because 
they engage electrons from vastly distant bands (similarly as the electron - electron inter- 
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action U) . Moreover, we have only the normal ordered part of this interaction entering to 
the Hamiltonian, and that should be small. The mean field value of this type interactions 
were already included by us in the fiow equation for /k,q (pT]). 

Having interactions shown in (p8D, one can induce from [t](1),5Hq (/)] next generation 
of higher order interactions, like for example c\^C2iCsiC4-^a}^aQ. They can be found in the 
standard canonical transformation as well (see H'^ in the Ref. |^). However let us repeat 
again, that in our case we get only the normal ordered forms of all such higher order inter- 
actions. So, hopefully, their infiuence on a physics of the effective Hamiltonian should be 
relatively negligible. 

V. CONCLUSIONS 



In summary, we formulated continuous canonical transformation for the double exchange 
model, eliminating from the initial Hamiltonian the exchange interaction term responsible 
for a violation of magnon number. Thus, in a limit / ^ oo, true magnons are obtained. 
Parameters of the effective Hamiltonian are determined via the set of the flow equations 
( p!^ - |T5|) and (^). Structure of the resulting Hamiltonian (P) is simple but the effective 



model parameters are computed selfconsistently taking into account effect of the higher 
order interactions (what is not possible in the standard single step transformation). These 
feedback effects are discussed by us on example of the magnon dispersion. We flnd a partial 
softening of the spin stiffness in a whole regime of hole concentration. 

Further studies are needed to solve the flow equations for a realistic 3D version of the 
DEX model. Another important issue not studied here is to extend the procedure to capture 
the damping effects for magnons and electrons. As discussed earlier in the context of different 
models Jll],^ one should study then the flow not only of the whole Hamiltonian H{1) but 
also for the particular operators like a^^l). Such an investigation for the DEX model is in 
progress and the results will be reported elsewhere. 
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